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Abstract 

We consider the motion of a particle in a force field subjected to adiabatic, 
fluctuations of external origin. We do not put the restriction on the type of 
stochastic process that the noise is Gaussian. Based on a method developed 
earlier by us [ J. Phys. A 31 (1998) 3937, 7301] we have derived the equation 
of motion for probability distribution function for the particle on a coarse- 
grained timescale At assuming that it satisfies the separation of timescales; 

<^ At <^ Tc, where Tc is the correlation time of fluctuations. \fJ-\^^ 
refers to the inverse of the damping rate (or, the largest of the eigenvalues 
of the unperturbed system) and sets the shortest timescale in the dynamics 
in contrast to the conventional theory of fast fluctuations. The equation 
includes a third order noise term. We solve the equation for a Kramers' type 
potential and show that although the system is thermodynamically open, 
appropriate boundary conditions allow the distinct steady states. Based on 
the exact solution of the third order equation for the linearized potential and 
the condition for attainment of the steady states we calculate the adiabatic 
noise-induced rate of escape of a particle confined in a well. A typical variation 
of the escape rate as a function of dissipation which is reminiscent of Kramers' 
turn-over problem, has been demonstrated. 
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I. INTRODUCTION 



About a century ago Einstein formulated the problem of Brownian motion as what is 
known today as stochastic processes. To explainEHi the motion of a Brownian particle on 
an observable macroscopic scale he introduced the coarse-graining of space and time in the 
dynamics. This has served as the basis for generation of successive levels of description from 
the microscopic to macroscopic realm in subsequent many other formulations. Although 
there exists no general program of coarse-graining it has been possible to describe the dy- 
namics fairly realistically in terms of some systematic separation of timescales consistent 
with experiments. 

In dealing with the stochastic processes one essentially examines the average motion of a 
system subjected to fluctuations which may be fast or slow depending on the correlation time 
Tc of the fluctuations compared to coarse-grained time scale, At over which one observes the 
average motion of the particle. While the stochastic processes with short correlation time are 
well understoodiHll, significant progress has been made in recent times in theories dealing 
with arbitrary long correlation time. However, one important constraint in this contextB 
needs to be emphasized. This is that the separation of timescales is just not enough and one 
has to specify further the nature of stochastic process by assuming it to be either Markovian 
or Gaussian or botfM. In fact these lie at the heart of the overwhelming majority of the 
traditional theories so far. The Ornstein-Uhlenbecki0'0 process is an age old standard 
paradigm in this respect. 

Very recently we have proposed a methodli'Ei for analyzing the stochastic dynamics 
of a particle subjected to external, adiabatically slow fluctuations where we do not limit 
ourselves to Markovian or Gaussian processes. It is apparent that to go beyond Markovian 
and Gaussian approximations the standard procedure of expansion of Master equationB 
(based on Markov approximation) leading to Kramers- Moyal [KM] equation or characteristic 
function methodElil (based on the calculation of moments and, in general, used for Gaussian 
processes) are not suitable for the purpose. Based on adiabatic following approximation we 
have derived the appropriate equations of motion for the probability distribution function 
on a coarse-grained timescale At assuming that At satisfles 

^ < At < Tc (1.1) 

where refers to the inverse of the damping rate (or the inverse of the largest eigenvalue 
of the 'unperturbed' system. This is in contrast to fast fluctuations characterized by the 
inequality Tc <^ At <^ -j^). therefore refers to the shortest timescale in the dynamics 
in the present investigation compared to short correlation time Tc of the fast stochastic 
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processes. The separation of the timescales (|1 . 1|) has been critically analyzed by Masoliver 
et. al.0 in their treatment of generalized Langevin equation for studying diffusion with 
external noise. 

In this paper we extend our earlier treatmentiilli to analyze and present a solution for 
the problem of motion of a particle in a force field subjected to external, adiabatic noise 
pertaining to the separation of timescales (|1 . 1| ) . The potential we consider here is of Kramers' 
type as shown in Fig.([l|). The physical motivation behind the solution is two fold : 

First, since it is subjected to fluctuations of external origin, the system is thermody- 
namically open (the close system on the other hand is characterized by internal noise which 
satisfies the fluctuation-dissipation relation). We search for the condition of an appropriate 
steady state for this open system. Similar studies for Gaussian processes have been carried 
out by Masoliver and othersi^'H. The study of open systems are specifically relevant for de- 
scribing the effect of pump fluctuations on the emission of a dye laseril, effect of fluctuating 
rate constants on a chemical reactionSil and effect of noise on parametric oscillatoril, etc. 

Second, once the condition for attainment of the steady state is realized, it becomes 
possible to consider the situation such that the particle in a force fleld, i.e., originally con- 
flned in a potential well may escape under the influence of external adiabatic noise by 
maintaining a steady state probability current over the barrier. It is therefore pertinent to 
calculate the rate of escape induced by this nonthermal activation in the spirit of Kramers 
and Smoluchowskiizhil and to elucidate the aspects of dependence of escape rate on dissipa- 
tion. The counterpart of the latter issue in the theory of fast fluctuations is the well known 
turn-over problemS^l. 

We thus intend to touch the three issues in the theory of adiabatic fluctuations 
in the subsequent sections : First, based on the method of 'adiabatic following 
approximat ions' together with a systematic separation of timescales (|1.1|) to carry 

out an expansion in as developed by us in two earlier paper JlUli we obtain an equation 
of motion for probability distribution function for a particle in a force fleld simultaneously 
subjected to external adiabatic noise. The correlations of fluctuating forces give rise to sec- 
ond and third order diffusion coefficientsEl. We discuss this issue in Sec. II. Second, although 
the system is thermodynamically open we look for the physically allowed steady states and 
show that application of the appropriate boundary conditions leads to distinct steady states. 
Sec. Ill is devoted to this aspect. Third, based on the exact solution of third order equation 
for the linearized potential and the condition for attainment of the steady states, we address 
the problem of escape of a particle conflned in a well in the spirit of Kramers-Smoluchowski 
theory. This forms the subject matter of Sec. IV. We conclude this paper by summarizing 
the main results and their experimental relevance in Conclusion. 
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II. MOTION OF A PARTICLE IN A FORCE FIELD IN PRESENCE OF 

ADIABATIC FLUCTUATIONS 



The equation of motion of a particle of unit mass in an one-dimensional extension where 
it is acted upon by an external field of force corresponding to a potential V{x) and an 
external, adiahatic stochastic force ^(t), can be written as follows 



x = ~V\x) + ^ai{t) . 



(2.1) 



Here we have considered the overdamped limit. (5 is the dissipation constant and a which 
is necessary to keep track of the order of perturbation, is a parameter determining the size 
of fluctuations of the external noise ^{t). We emphasize here the two points : (i) we do 
not put restriction on the type of stochastic process ^(t) that the noise is Gaussian. This 
has attracted so much attention in the recent literature that it is necessary to point out 
that no such assumption has been made. The only restriction we make on the nature of the 
stochastic process ^{t) is that its correlation time Tc is very long, i.e., it corresponds to the 
separation of the timescales implied in the inequality (|1 . 1| ) . Also note that the inequality 
implies the overdamped limit, (ii) We assume for convenience, without any loss of generality, 
that {iit)) = oi. 

In a preceding paperEl we have derived the equation of motion for probability density 
distribution function P{x, t) in phase space corresponding to the Langevin description (2) 
where the associated timescale satisfies the inequality ( |1.1| ). We have shown that P{x,t) 
obeys the differential equation of motion which contains third order terms (beyond the 
usual Fokker-Planck terms) leading to non-Gaussian noise. The appearance of these terms 
is generic for the stochastic process we consider here. The general expression for time 
evolution of probability density function is given byS 



dP{x,t) 



dx 



/■oo 

V-Fo + a^ V- / {FiV.r-Fi{x'^)) 
Jo 

r (FiV_.-Fi(x-0) 

WV- r {FiV^r ■ F, 
Jo 



dx 



dx 



dr 



dr 



dx 



dx ^ 




dx 





P{x,t) 



(2.2) 



where Fq and Fi refer to the unperturbed and the fluctuating terms, respectively, corre- 
sponding to Eq.(2) as given by 



^0 = -^y'i^) 



and Fi 



am 



(2.3) 
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The symbol V is used for the operator that differentiates everything that comes after it 
with respect to x. V_r denotes the differentiation with respect to X-r- The Jacobian 
defines the mapping x for the unperturbed motion and is given byEl 



dx 



dx 



dx 



dx 



1 + -V"{x)r + 0{T^) 



(2.4) 



where it has been assumed that x varies very httle in r (in the scale of Also explicitly 
we have 



d_ 

dx ' 



(2.5) 



and 



F„(x-^) = -jV'(x-) = -i [V\x) - T V"[x)\ 



(2.6) 



(2.7) 



Making use of the relations (4-8) and after collecting the terms of the order of a^, all the 
four terms in Eq.(3) can be simplified further to obtain [ some details are outlined in the 
Appendix- A ] : 

1 



dP(x,t) 1 d ,,,,, 2 1 d'^P(x,t) 2 

Co, Ci and C2 in Eq.(9) are given by 



dx^ 



[V'{x)P{x,t)] . (2.8) 



Coi 
Co 
Cl 



Co - Cl 

oo 

{^{t)at-r))dT 







Co 



Tcieit)) 

POO 

C2= / Tmat~r))dT . 

Jo 



(2.9) 



We now put a = 1 for the rest of the treatment. 

The above equation describes the time evolution of an overdamped particle in a force 
field (derivable from a potential V{x)) simultaneously subjected to an external adiabatic 
stochastic force. Cq, Ci and C2 measure the strength of the noise term. While the first term 
in Eq.(9) can be identified as the usual deterministic dynamical term, the second and the 
third terms refer to second and third order diffusion coefficients due to stochasticity C,{t). 
The remarkable departure from the standard form of Fokker-Planck equation (Smoluchowski 
equation) is due to the presence of the third order noise. 
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Let us now digress a little bit about Eg. ( p.8|) which forms the basis of this paper. The 
very appearance of third order term apparently suggests an immediate bearing with the 
familiar KM expansion. We emphasize that Eq. (|2.8|) is not a truncated KM expansion. The 
KM expansion is based on an expansion of Master equationEl in a series in r and the terms of 
(9(r^) are neglected. Or in other words the coefficients in a KM expansion originate from the 
corresponding moments ( of the distribution function ) which are assumed to be linear in r. 
This severely restricts the use of higher order moments for finite r. This is mathematically 
formalized in terms of Pawula TheoremEl ( whose proof concerns the relation between the 
moments pertaining to the KM expansion only. Note that the familiar Wigne equation 
for probability distribution for a cubic potential is a third order equation which has nothing 
to do with a KM expansion or Pawula theorem ) which precludes the occurrence of nonzero 
terms beyond Fokker-Planck in a KM expansion. Thus if there is any finite third order term 
it has to be treated as a perturbation. In many situations one encounters serious interpretive 
difficulties since the probability distribution functions often turn out to be negative. 

The expansion leading to Eq. ( |2.8| ) on the other hand is based on an expansion [ Eq.(9) 
oHl in af3 ^ where a is the strength of noise and f3 is the damping constant. Note that f3 ^ 
defines the shortest time scale in the dynamics. That the expansion is in a(3~^ is evident 
from the appearance of ^, ^ and -p- factors in the successive terms in Eq.(p.8|). Although 
the convergence of formal KM expansion, in general, is not guaranteed ( it is built in from 
outside as done by van Kampenl in terms of expansion ) it has been proved thatEl our 
expansion scheme is convergent in the adiabatic following limit. The appearance of third 
derivative term is generic since it appears in the same order ( a^-order ) as the second 
derivative term in the present expansion schemed and is characteristic of non-Gaussian 
features. Thus the third order term cannot be treated as a perturbation. 

Eq. (|2.8| ) is derivedlil'lli on the basis of two approximations : (i) adiabatic following ap- 
proximation which results in a convergent expansion in ( valid for very slow noise 
processes corresponding to the timescales (|1 . 1| ) ) which is complementary to the cummulant 
expansion based on an expansion in ar^ ( valid for very fast noise processes! ). (ii) decoupling 
approximation. We have calculated the error caused by decoupling and shownEl'S that it 
is of the order of a'^(3~^. The corresponding error in decoupling in the case of fast fluctu- 
ations is a'^Tc. We explicitly point out that no other approximation is required to develop 
the theory further as done in this paper and given the appropriate boundary conditions the 
derived solutions are well behaved positive definite probability distribution functions. 
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III. THE SOLUTION OF THIRD ORDER EQUATION FOR A KRAMERS' 
TYPE POTENTIAL : STEADY STATE PROBABILITY DENSITY 



A. The general solution 

To start with we now recast the third order equation (9) in the form of the famihar 
continuity equation and identify the current S{x,t) as follows; 



dP{x,t) _ dS{x,t) 
dt dx 

where 



(3^1) 



Si..t)^-lv'i.)P-f'^^f,^lV'i.)P} . (3.2) 



In the steady state one puts 

dP{x,t) 



(3.3) 



dt 

to obtain the following steady state current J as given by 

./^-iv'WP.-f f + . (3.4) 

where Pg is the steady state probability distribution. Henceforth for the sake of brevity we 
omit the subscript s from for all the subsequent calculations to follow and denote P{x) 
as the steady state probability distribution function. 
Multiplying both sides of Eq.(14) by P^/c2 we obtain 

^iV'(.)P}-JJ^-HV'MP} = ^-^ . (3.5) 
where a and b are given by 

a = and b = — . (3.6) 

We now explicitly make use of the Kramers' type potential V{x) as shown in Fig.(||) for the 
problem. We consider by linearizing it at x = 0, 

V{x) = Eo- \ul x" (3.7) 

where ujq refers to frequency of the inverted well. Eq defines the potential at the barrier top. 
The probability distribution function P[x) therefore satisfies the following equation 
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(pp 



'2ujI + a\ dP , „ 
° ' -hxP 



dx 



C2 <^o 



Defining 
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and D 



C2 ci;5 



the Eq.(18) reads as 



d^P dP , ^ ^ 

x-rr + 1 b X P = D 

dx'^ dx 

It is convenient to make the following substitution 
in Eq.(20) to obtain 



X 



dx"^ 



+ x 



dW 
dx 



i^-lf + hx' 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



Let 



We then rewrite Eq.(22) as follows, 

x^-— + X — (u^ + bx^)W = D x^^" 

dx^ dx 

Prom the definitions of 7 and a (Eqs.(19) and (23)) we have 



2 V c^gcj 



(3.13) 



(3.14) 



(3.15) 



The structure of Coi(= cq — Ci) suggests (see the definition (10)) that by virtue of adiabatic 
stochasticity, Ci is much smaller compared to Cq and Cqi is always positive. This ensures that 
V as defined in Eq.(25) is always positive. 

It is convenient to make a further substitution for independent variable x as 



X 



This reduces Eq.(24) to 

,d?w 



+c 



dW 



(i.2 + (2) ^ 



D 



(3.16) 



(3.17) 
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The homogenous counterpart corresponding to the above Eq.( p.l7| ) is the standard modified 
Bessel equation of order v. The general solution of Eq. ( p.l7| ) can be written as 



W{0=AI,{C) + BK,{C) 



where luiC) and K,y{() are modified Bessel functions of order z/; A and B are the two 
arbitrary constants of integration corresponding to the homogenous part of Eq.f p.lTl) . The 
D containing term results from the particular integral corresponding to the inhomogenous 
contribution of Eq.( |3.17D obtained by the method of variation of parameters. Making use 
of the relations (26) and (21) we revert back to the original variables x and P{x) to obtain 
the general solution of Eq.(20) as 



5X 



P{x) = A x-" luiVbx) + B x-"" K^iVb: 

I^iVbx) f K^{Vbx') dx' - K^iVbx) f h 



+D x"" 



. (3.19) 



B. The boundary conditions and the normalized probabihty distribution 

At this juncture it is necessary to specify the boundary conditions. We impose the 
following natural boundary conditions on the solution (29) 

(i) P{x) vanishes for |x| — > oo 

(ii) P{x) remains finite at a; = 0. 
To this end we proceed as follows : 

(i) First we note that a modified Bessel equation of the form 

x^ y" + xy' - {x^ + y = (3.20) 
has an irregular singular point as x — oo. The leading behavior of the solutions areS 

y{x) = h{x) ~ Ci x-^l'^ , x-^oo (3.21a) 

and 

y{x) = K^{x) ~ C2 x-^''^ e""^ , a; ^ 00 . (3.21b) 

We thus observe that Iy{x) diverges exponentially for large x. By applying the boundary 
condition (i), i.e., P{x) vanishes for large x we see that the constant A in the general solution 
(29) must be zero. Therefore we have 
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P{x) = B x^" KiVbx) + D x'" I^iVbx) J K^{Vbx') dx' 

-D K^iVbx) j^^^'" h{Vhx') dx' . (3.22) 
B and D are the two remaining constants to be determined. 

(ii) We now turn to the second boundary condition, i.e., the finiteness of P{x) at x = 0. For 
a fixed (z/ > 0) we know thatS 

K^{x) ~ ]-T{v) (]- x^ , for a; ~ (3.23) 



2 ' ' V2 

which imphes that Ky{x) has a singularity at x = 0. This reveals that P{x) is singular at 
X = 0. Our strategy here is to remove this singularity by having an appropriate relation 
between the two remaining constants B and D in the solution (36). To derive this relation 
we proceed as follows: 

We first derive the asymptotic expansion for the solution of the modified Bessel equa- 
tion (30) whose leading behavior is (31a). To do this we peel off the leading behavior by 
substituting 

y{x) = Ci x-^/2 gx ^^^^ ^3 24) 

into the modified Bessel equation (30). The equation satisfied by w{x) is 

x^ w"{x) + 2x^ w'{x) + - w{x) = . (3.25) 

We seek a solution of this equation of the form w{x) = 1 + e{x) with e{x) ^ 1 (x — oo). 
e{x) satisfies the equation 

x^ e"{x) + 2x2 e'{x) + - '^^^ + " '^^) = ' (3-26) 
which may be simplified by the approximations 

Q - z/2^ e < i ~ ' ^ ' ^ ^ °o • (3-27) 

We make the second of these approximations because we anticipate that e decays like a 
power of X as X — > oo. The resulting asymptotic differential equation is 

2x2 e' ~ (u'^ — t] , X ^ oo 



Ordinarily the solution to this equation would be e{x) ~ c as x — > oo where c is an 
integration constant. However since e{x) ^ 1 as x — > oo we must set c = 0. The leading 
behavior of e{x) is then given by 
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1 1 2\ -1 

B{x} ^ \ - — -y ] X ,x— i>oo. 



This kind of analysisE3 can be repeated to obtain all the terms in the asymptotic expansion 
of w{x) as X — s> oo. However the leading behavior of e{x) suggests that w{x) has a series 
expansion in inverse power of x. Thus to simplify the analysis we assume at the outset that 



oo 
n=0 



(x — > OO , ao = 1) 



(3.28) 



Substituting this expression into the differential equation for w{x) gives 



n{n + 1) a„ x " - 2 ^ n a„ " + ( ^ " ) X! ^" ^ " ~ 



n=0 



n=0 



(3.29) 



n=0 



Since the coefficients of any asymptotic power series are unique we equate to zero the coef- 
ficients of all power of - in the above relation 



X 



n-\ — 
2 



a„ — 2 (n + 1) a„+i = , n = 0, 1, 2 . . . 



(3.30) 



Solving this recursion relation and using ao = 1 we obtain 

(41/2-12) (4z/2-l2)(4z/2-32) 

w{x) ^ 1 - ^ ^. „ - + — ' 



1! 8x ' 2! (8x)2 

From the ratio test we see that the radius of convergence R of (41) is 

2(n+ 1) 



X ^ oo 



(3.31) 



R = lim 



lim , ,„ „ 

n^oo _|_ ly — 







(3.32) 



unless the series (41) terminates, which it does when 

1 3 5 

z/ = ±- , ±- , ±- , . . . . 

2 ' 2 ' 2 ' 

When this happens, the finite series (41) when multiplied by e~^/ y/x gives an exact solution 
to the modified Bessel equation. 

Similarly, the complete asymptotic series for the function whose leading behavior is given 
by (31b) is 



y 



(x) ~ C2 x-^/2 e"" w{x) 



where 



w[x] 



^ , (4Z/2-12) ^ (4z.2-l2)(4,.2_32) ^ 



1! 8x 



2! (8x)2 



X — * 00 



(3.33) 
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By global analysisE3 the two constants Ci and C2 can be derived as 



Ci = (27r) 



"1/2 



C2 = \W2 , (3.34) 

to get the modified Bessel functions /^(x) and Ky{x) respectively. 

Thus, in general, /^(x) and K^{x) are represented by infinite power series. The above 
analysis reveals that for u = ±1/2, ±3/2, . . ., the series terminates and from the asymptotic 
expansion we get an exact solution. In such e., when the power series terminates it is 

possible to find a relation between the two constants B and D of the general solution ( p. 22] ), 
such that the singularity of P{x) at a; = is removed. Thus both the boundary conditions (i) 
and (a) [i.e., P{x) = for |x| 00 and P{x) is finite at a; = 0] are satisfied by the solution 
( 15*. 22[ ) provided u is an odd half integers, i.e., 



1 

z/ = n + - , n = 0,1,2,... . (3.35) 



By Eq.(25) we have 



1 

'^+2 ' ^ = 0' 1'2, 



But n = implies cqi = or cq = ci which is not allowed physically. Therefore we have 

£^ = n; n = 1,2,3,... , (3.36) 

zuJq C2 

In passing, we mention here that the integers n characterize the distinct physically allowed 
steady states of the thermodynamically open systems. 

With the above mentioned restricted values of u = n + ^ , n = 1, 2, 3, . . .) we now 
explicitly calculate the several quantities [ details are given in Appendix-B ]which appeared 
in the general solution P{x) in Eq.( p.22[) . We finally obtain : 



Pn+^i^) ~ ^ Y 2 51/2 ^ fk ^k+n+l 



k=0 

-5^ EE g {(-)■ + (-1)-} /r f: 4^ (3.37, 

where we have defined 

^ [n + ky. 

2kh^/^k\ {n-k)\ ■ ^ > 

We are now in a position to remove the singularity of P(x) at x = . To achieve this we 
seek a relation between the constants B and D . For this we now expand the exponential 
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Eq.(47) . It is convenient to write a few nontrivial steps explicitly after this expansion. 
3 thus have 



m 

We thus have 

n oo 



fe=0 a=U 

71 " " "^'^ (n — kW 1 1 

7-) n n n-k (n — hW 1 1 

^ V f-lV'"" V ^ 



/c=0 a=U 

D " """^ fr? - 2V 1 1 

n " " r?! 1 1 

D " fn - IV 1 1 

(where we have carried out the summation over k) 



D " 



" , n' 1 

E(-i)' /o -r 

j=0 



n-l 



+ E(-i)'/fc"/r 



(n-l)! 1 1 



3=0 



h 



" r?' 1 1 

E(-ir"A"/o -7 



Ct? - IV 1 1 



(3.39) 



(where we have replaced the dummy index i by A; in the D-containing sums) 

To remove the singularity at x = from the above expression for P„_,_i(x), B must be 
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related to D in such a way that all powers of ^ vanishes identically. For this the coefficients 
of -jTFT from the summation are equated to zero. This gives 



Note that the relevant coefficients contribute from the 5-containing term for a = 1 and from 
D-containing term for j = n. Putting the value of from Eq.(48) in Eq.(50) we obtain the 
desired relation between B and D. 

r;' 1 

Prom Eq.(51) it reveals that the relation between B and D is independent of any dummy 
index i, j or k and only depends on n. Thus for any n the relation is unique. With this 
choice the probability P„+i(x) (Eq.(47)) becomes finite for all x and is given by 

-in -1 



PnHi^) - (-1)" D 2- n\ E e-^"'^^. 



k=0 



EE E{(-1)^ + (-irn /r ^^V^ 7^^:3-717 (3-42) 



The above probability distribution function which remains undetermined upto the constant 
D (which we are to determine shortly) vanishes as x — > oo and remains finite at the barrier 
top X — 0. For illustration, we explicitly write some of the probability distribution functions. 
For n — 1 case we have, 

P,.(.) ^ -2D r=/= e-V-- (i, , ^) _ f (i (3.43a) 

X - (3.43c) 



where the probability distributions for x < are obtained by noting the symmetry of the 
differential equation. 

To normalize the probability and thereby to calculate the current over the top of the 
barrier located at a; = 0, we consider the condition for normalization: 

/+00 
P,,^{x)dx^l . (3.44) 
-co 

Since the probability decreases rapidly (due to the presence of the modified Bessel function 
K^{x)) we approximate the integral (54) by 
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r+A 

/ (a;) dx = 1 ] A large but finite (3.45) 

J- A 2 

where ±A approximately indicates the two zeroes of the inverted potential V{x). By Eq.(17), 
A may be expressed in terms of Eq the height of the barrier as, 

A = (2Eo/^o')'/' • 

By symmetry we rewrite 

.A 

/ P^+iix) dx = 1/2 ; A large but finite . (3.46) 
Jo 2 



Due to the presence of the powers of ^ in the expression for P^^i{x), it is not possible 
to carry out the integration in Eq.(56) in a straightforward manner. We therefore employ 
a limiting procedure and make use of the expression (52) for P„_,.i(x) in (56) to obtain an 
expression for the following normalization constant D 

D = (3 47) 

2 Ein(v^A) + 2^^+! n\ ELoE^io"^ 2^ iJ-~,Zn^,-,) ' 

where A is large but finite and the function Ein(x) is defined a; 

Ein(x) = £(-1)^ 1^ . (3.48) 

The expression (52) together with the normalization constant D as given by (57) 
yields the complete and exact analytical expression for the probability distribution func- 
tion Pn^i{x). In Fig.(0) we draw two typical normalized probabilities for n = 1 and n = 2 
for the parametric values (3 = 1.0, A = 2.5 and C2 = 0.9. The effect of third order noise 
is illustrated in Fig.(|^), where we exhibit the probability distribution functions for several 
values of the third order noise strength C2. One observes that the distribution gets flattened 
as C2 increases as expected. 



IV. DYNAMICS OF BARRIER CROSSING INDUCED BY ADIABATIC NOISE 

We now extend the above analysis to elucidate the following problem of dynamics of 
barrier crossing. 

An overdamped particle moves in an external fleld of force and in addition to this is 
subject to an adiabatically fluctuating force of external origin. The conditions are such that 
the particle is originally caught in the potential well in V{x), but may escape in course of 
time over the barrier. Our object here is to calculate the rate of escape from this well. The 
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analogous problem for the case of fast fluctuations is the celebrated Kramers' problem of 
Brownian motion in phase spaced. Our calculation rests on the third order equation of 
motion (9) obeyed by probability distribution function derived in Sec. II. 

To proceed further we employ the The popular fiux-over-population method originated 
by Farkasii many years ago. The calculation rests on the evaluation of two quantities; (i) 
the steady state current J over the barrier top, (located at x = 0) that results if the particles 
are continuously fed into the domain of attraction (say, in the region of left well) and are 
subsequently and continuously removed in the neighboring domain of attraction, (ii) Steady 
state population Ua in the initial domain of attraction, i.e., the left well. The rate is defined 
by 

K, = J/na . (4.1) 

The method has been used by Kramers in his seminal work on barrier crossing dynamics 
and many others over the several decade&i. 



A. Calculation of 

For calculation of it is necessary to evaluate the stationary probability density Ph{x) 
near the bottom of the well corresponding to a zero current (J = 0) situation along the x 
co-ordinate. 

We first linearize the potential V{x) as shown in Fig.(|I]) near the bottom of the left well, 
at a; = —A, so that we approximate 

V{x) ^ \ ul (x + Af , (4.2) 

where Ub refers to the frequency at the bottom of the left well. Eq.(60) and J = condition 
reduce the third order equation of motion (Eq.(9)) for probability density Pb{x) inside the 
well to the following form: 

{x + A)^ + i^~b{x + A)Pb = . (4.3) 
The above equation is valid near the bottom of the left well (x ~ —A). Here 7' is defined as 

y = ^. (4.4) 

where a is as given by Eq.(16). Putting z = x + A and y = zPb, Eq.(61) can be transformed 
as follows : 

z'^-{2-^')zf^ + [i2-y)-bz']y = . (4.5) 
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From Eqs.(62) and (16) we have 

2-7' = ^. (4.6) 
Note that 2 — 7' is a positive quantity. We write 

2 -7' = a . (4.7) 

Therefore Eq.(63) is given by 

Substitution of y{z) = z^^""^^^ W{z) in Eq.(66) yields 
where 

^' = ^ • (4.10) 

The solutions of Eq.(67) are again the modified Bessel functions, I,y'{Vbz) and Kj,t{Vbz). 

Reverting back to original variables, the general solution for the steady state probability 
distribution near the bottom of the left well is given by 

Pbix) = A'{x + Ay' h,[Vh (x + A)] +B' {x + AY' K,\Vh {x + A)] , (4.11) 

A' and B' are the two arbitrary constants of integration. 

By demanding that Pb{x) must vanish at infinity, we require 

A' = (4.12) 

and therefore 

Pi,{x) = B' {x + AY' K,.[Vh {x + A)] . (4.13) 

Furthermore, we note that although K^i itself is singular at x = —A the presence of {x + AY' 

in Eq.(71) assures that the probability (71) remains finite at x = —A. To verify this 

assertion, we use the property of modified Bessel function K^i{z) for fixed v' and 2; — 0. 
Kyi{z) behaves asii 

K^,{z) ^\t{v') , Rez/' > . (4.14) 



2 ' ' V2 
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Hence 

Pb(-A) = B' lim {x + Af K^,[Vb{x + A)] 

= B' ^^i^^ . (4.15) 

We thus see that the steady state probabihty P\y ( given by Eq.(71) ) is finite at the bottom 
of the left well. 

The above solution Ph{x) must now be subject to the following boundary condition : 

A(-A)=Pr^(-A) (4.16) 

n+- 

where the stationary probability '^{—IS) corresponds to the vanishing current J = 
along X pertaining to the homogenous version of Eq.(20). As usual, Pi ^ (x) must also 
satisfy the boundary condition that for |x| oo, ^(x) vanishes. Such a solution is 
immediately apparent from our earlier analysis of Sec. III. Thus 

Pt'-Hx) = B E A" ; ^ > (4.17a) 

where is as defined in Eq.(48). 

Making use of Eqs.(73) and (75) in Eq.(74) we obtain a relation between B and B' . 

^' = ^ V 26V^ T{v') 2-'-i ■ ^^-^^^ 

Therefore P6(x) in Eq.(71) may be expressed as 

^^(^) = B ;n ___(, + A).' ;^^,[V6(x + A)] . (4.19) 

The above distribution which is valid near the bottom of the left well may be used to 
calculate the population inside the left well as, 

na = 2 f Pb{x) dx . (4.20) 

Due to the presence of Ky{x) the probability Pb{x) is a rapidly decreasing function. We may 
extend the above integration limit to infinity. This yields ( using Eq.(71) )lli 

Ua = B' ^ ^ ^ . 4.21 

Using the relations (76) and (79) we finally have 

na = V2.B A_±l2 Jlyir--' ^ r3/^ . (4.22) 
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B. Calculation of escape rate 



Having determined the population of the left well and the steady state current J from 
Eqs.(80) and (19), respectively we are now in a position to calculate the escape rate in terms 
of the Eq.(59). We thus obtain (for calculation of J the linearization of the potential V{x) 
at X = which results in Eq.(17) has been carried out) 



C2UjI (D\ r(z/') J 6^/4 



(4.23) 



^n., = i-ir' 4^ ^ { , : ^ \ ■ (4.24) 



2 



Making use of the relation (51) between the two constants B and D in Eq.(81) one finds 

The expression (82) can be simplified further by noting the following relations. First, we 
have from Eqs.(19) and (23) 

1/ = i (7 - 1) where 7 = 2 + -^ , (4.25a) 
2 

and from Eqs.(62), (65) and (68) 

z/' = - (1 - 7') where 7' = 2 - 4 • (4-25b) 
2 uJh 



Thus we have 

2^ I ujlul 



v + v' = t.a (^^V^ . (4.26) 



Furthermore Eqs.(16) and Eq. ( |3.36| ) suggest that 



C2 

Use of Eq.(85) in Eq.(83b) yields 



a = 2nujl ; h=— ; n = l,2,... (4.27) 



^2 ^ 

z/' = n-|-- ; n = l,2,... (4.28) 



The above relation together with the expression for (Eq.(48)) leads to the formula for 



the transition rate /C„_,_i as follows 



,2 



„ , ,2 1 T{n'^-\) { p-VbA 



19 



A-s refer to the zero's of the potential V{x) as shown in Fig.([^) and by virtue of hnearization 
of V{x) at X = (Eq.(17)) A is approximately given by 



The above expression can be made more transparent by demonstrating a representative 
transition rate, say, for n = 1 as follows : 



The above expressions are analogous to Kramers' formula for the rate of escape from a 
potential well over a finite barrier of height Eq under the influence of an external nonthermal 
adiabatic noise. What we have shown here is that within a linearized description of the 
potential, the corresponding diffusion process can be described exactly for dissipation (3 
pertaining to the timescale ^ <^ At <^ Tc . 

The escape rate expressions derived above suggest that the rate approaches zero both for 
(3 oo and /5 — > 0. This behavior is somewhat reminiscent of Kramers' theory, where it was 
noted earlier that these two limiting behaviors imply a maximal rate at some damping value 
f3. The rate therefore undergoes a turnover in a form of a bell-shaped curve. In Fig.(|^) we 
plot a representative variation of the escape rate versus dissipation f3 for different third order 
noise strength. With increasing friction, the rate undergoes a turnover from an increasing 
behavior at low friction to an inverse behavior in the high friction limit. 

Since the driving noise is of nonthermal origin, the escape rate expression is devoid 
of any temperature. Temperature is characteristic of a closed thermodynamic system at 
equilibrium. What we have here instead is the ratio of strength of nonthermal noise to 
dissipation (3 in the exponential factor of the rate expressions. As mM pointed out that in 
the steady state such a parameter might play the role of temperature in the open system 
which characterizes the steady state. 

As emphasized earlier that since the noise is external and the noise and dissipation have 
no common mechanistic origin (in contrast to what one observes in theory of Brownian 
motion) the steady state is not allowed arbitrarily. The system can attain the steady states 
depending on the specific integers n {n = 1,2,...) which uniquely connect the parameters 
Co, ci, C2 according to the relation derived in Sec. III. 



In this paper we have presented a solution for the problem of motion of a particle in 
a force field simultaneously subject to an external adiabatic noise characterized by long 




(4.30) 




(4.31) 



V. CONCLUSIONS 
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correlation time without keeping any restriction on the type of noise that it is Gaussian. 
Specifically we have calculated the rate of escape of the particle over the barrier initially 
confined in a well induced by nonthermal fiuctuations. The theory rests on a perturbative 
expansion in \fi\~^ ( where \fi\ is the damping constant or the largest eigenvalue of the 
unperturbed system ) pertaining to the separation of timescales (|1.1|) as carried out in our 
earlier papergl3li2l. The main conclusions of our study follows: 

(i) An analogue of Smoluchowski equation in the case of adiabatic non-Gaussian noise 
processes (Eq.(9)) has been proposed. 

(ii) The undergoing stochastic process is characterized by third order noise which is 
responsible for non-Gaussian featureslil. 

(iii) Given the appropriate boundary conditions the third order equation admits of exact 
solution for the linearized potential. 

(iv) The interplay of the characteristic linear dissipation of the system and the external 
noise leads to physically allowed distinct steady states subject to appropriate boundary 
conditions. 

(v) In the spirit of Kramers theory we have solved the problem of escape of the particle 
confined in a well and have shown that the escape rate exhibits a turnover as one passes 
from the relatively low dissipative to the strong dissipative regime ( Eq.(87) ). 

In conclusion, we have thus discussed a number of basic issues in the classical theory 
of motion of a particle in a force field in presence of external, adiabatic fiuctuations. In 
view of several experimental investigations on external noise-induced processes in the past, 
the study of thermodynamically open systems has been specially relevant in various ar- 
eas of physical and chemical sciences. We particularly mention the following examples : 
A dye laser with fiuctuating pump parameter shows interesting qualitative changes in the 
stationary distributionsil. The shape of the distribution also changes as a function of in- 
creasing fiuctuation strength in the case of a two-species chemical reaction with a fiuctuating 
rate coefficienl0li. The early work on electronic parametric oscillator driven by external 
noise which exhibits the transition from non-oscillatory to oscillatory behavior is also worth- 
mentioningii. Although the driving noise processes in the above mentioned cases are fast, 
suitable extension to adiabatic noise limit (such a typical case had been discussed by us ear- 
lier in0 in detail in connection with population inversion in a two-level atom by adiabatically 
varying the field strengthil^ might lead to experimental situations which are relevant to 
the present theoretical context. We believe that studies in this direction are worth-pursuing. 
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APPENDIX A: SIMPLIFICATIONS OF THE TERMS IN EQ.(3) 



Making use of the relations (4-8) we simplify below the four terms as appeared in Eq.(3). 
First term : 



- V • Fo P{x, t) = -^[-^ V'{x) P{x, t)] 
= l^[n^)P(.,.)]. 



(Al) 



Second term : 



a 



V-^"|(FiV_,-Fi(0) 



dx 



dx 



dr} P{x,t) 



d 



dx Jo P 
1 92 



1 



1 - -V"{x-nr 



1 + i V'\x)r 



dr P{x,t) 



— a 



— a 



(3"^ dx^ Jo 
1 92 



/?2 dx'^ 



/ mat-r))[i-o{T')]dTP{x,t) 

Jo 

poo 

/ {mCit-r))dTPix,t) 
Jo 



(A2) 



Third term : 



r (FiV_.-Fi(x-^)) 



dx ^ 



dx 



o 1 d^P 



1 - -V"{x-^)t 



drj P{x,t) 

d 1 dm 



dx P dt 



(t-r) 



1 + i V"{x)t 



Fourth term : 



dr P{x, t) 

(A3) 
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dx 



1 - ^l^"(x-nr 



drj P(a;,t) 

1 - ^V"{x-nr 



d_ 

dx 



[V'{x) - tV"{x)]) 



1 53 ,00 



1 



1 + - V"{x)t 



dr P{x,t) 



—a 



1 



/?3 dx'^ Jo 
1 (9^ 



rmat-r)) [V'{x)-tV"{x)] [1 + -V"{x)t] dr P{x,t) 



(A4) 



APPENDIX B: CALCULATION OF VARIOUS QUANTITIES FOR 
SIMPLIFICATION OF EQ.(32) 

We explicitly calculate the several quantities which appeared in the expression for P{x) 
in Eq.( |3.22D . These are necessary for the derivation of Eq.(47). To this end we first note 
that 

T ( /h\- 1 ^ V- (-1)' {n + k)\ 
Using the expression for (Eq.(48)) 

(n + A;)! 



rn 
Ik 



2k 5fc/2 ^1 (^n-ky. 

we rewrite 



fc=0 



1 pV^fta: 

^-i(^-) = 7J7?7jS(-i)'A"^ . (B4) 



fc=0 ^ ^ 

Therefore we have 



= J—- Y: fk e-^^ X--' (B5) 



and 
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K^{Vbx) dx 



2 61/2 



fc=0 



Integrating successively by part^l (n — k) times we obtain 



/VOX 
x'' K^iVbx) dx = - 



2 61/2 



Vbx ^ ^ yra 



A;=0 i=0 



j! (v^)"-fc-i+i 



Similarly we have 



/VOX 
x'^ Iij{yhx) dx 



n 11— k 



fc=o i=o 



n fu {n-k)\ 



X-' 



V2 TT 61/2 

Hence from (B7) and (B4) we obtain 

x"'^ h{^x) J'^^" KiVbx) dx 



j! (v^) 



X 



2 61/2/ 



fc=0 j=0 



1 " pVfca:: 



1 n n n-k (n — hW ^j-i-n-l 

■ A E E E(-ir /r A" 

zv j—n k=n i=n J- 



j=0 fc=0 j=0 



(V6)" 



Similarly, 



1 



x'^ Iy{\/hx) dx 

n n n-k {n - k)\ x^-^"^-^ 



^EEE(-ir"/r/."— : 

^VO j=o fc=n 1=0 J 



n—k—j+l 



(B6) 



(B7) 



(B8) 



(B9) 



(BIO) 



i=0 k=0 j=0 

The expressions (B3,B4) and (B9,B10) can now be utilized in Eg . (13.221) to obtain the ex- 
pression for probability distribution function P(x) ( Eq.(47) ). 
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FIGURES 



FIG. 1. A schematic plot of the Kramers' type potential V{x). 

FIG. 2. The normalized probability distribution function -P„_|_i(a;) is plotted as a function of 
X forra = 1 and n = 2 ( /3 = 1.0, A = 2.5 and C2 = 0.9 ). 

FIG. 3. The normalized probability distribution function ^3/2(2^) is plotted as a function of x 
for various values of the third order noise strength C2 ( /3 = 1.0 and A = 2.5 ). 

FIG. 4. Escape rate /C3/2 is plotted as a function of the characteristic dissipation /3 of the 
system for various values of C2 ( cqi = 7.0, loi, = 0.80 and A = 2.5 ). 
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